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ABSTRACT 

A  time  domain  measurements ,  electromagnetic  scattering 
laboratory  has  been  developed  to  provide  a  means  of 
experimentally  measuring  the  transient,  scattered  electro- 
magnetic fields  from  symmetric  objects  for  impulse  response 
determination  or  target  classification.   The  evolution  of 
the  physical  laboratory  is  described  in  terms  of  the 
development  of  the  transmitting  and  receiving  antennas ,  the 
interface  of  the  data  processing  hardware,  the  signal  pro- 
cessing scheme,  the  software  routine  functions,  and  the 
fabrication  of  scattering  range  targets . 

The  smoothed  impulse  response  of  a  thin  wire  scatterer 
is  experimentally  determined,  and  a  comparison  of  these 
measurements  is  made  to  a  theoretical  calculation  of  the 
response . 
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I.   INTRODUCTION 

A.   BACKGROUND 

Since  the  early  nineteen  sixties,  research  in  time  domain 
or  transient  electromagnetics  has  continued  to  expand  into 
new  areas  of  application.   Current  research  involving  tran- 
sient electromagnetics  includes  wide-band  radar  design, 
aircraft  and  missile  radar  cross  section  reduction,  as  well 
as  the  evaluation  of  the  vulnerability  of  ships,  aircraft, 
and  ground  installations  to  high  level  nuclear  electromag- 
netic pulse  (EMP)  phenonmena.   A  direct  time  domain  transient 
solution  to  an  electromagnetic  interaction  problem  can  offer 
physical  insight  not  readily  apparent  from  a  solution  ob- 
tained by  classical  frequency  domain  techniques.   Reasons 
to  support  this  are  that  the  scattering  impulse  response  of 
an  object  is  related  to  its  geometrical  configuration;  the 
radar  cross  section  can  be  obtained  from  the  impulse  response 
by  a  Fourier  Transform,  and  the  scattering  response  of  an 
object  to  any  waveform  may  be  obtained  from  the  impulse 
response  by  convolution. 

As  a  practical  application,  advances  in  radar  target 
identification  brought  forth  by  using  transient  response 
techniques  could  be  profound,  indeed.   Consider,  for  example, 
the  military  and  commercial  impact  if  a  radar  operator  could 
positively  identify  a  radar  return  on  his  indicator  as  an 


F-M-  or  MIG-21  aircraft  and,  in  addition,  also  be  able  to 
determine  the  load  configuration  such  as  missiles,  bombs, 
external  fuel  tanks,  or  electronic  countermeasures  gear. 
This  idea  is  conceivable  with  advances  in  wide-band  radar 
technology  based  on  transient  electromagnetics  principles. 
Generating  an  electronic  signal  which  approximates  an 
impulse  was  not  possible  until  recent  developments  in 
subnanosecond  pulse  technology.   This,  coupled  with  the 
advances  achieved  in  high  speed  digital  processing,  gave 
impetus  to  the  evolution  of  the  methodology  necessary  to 
solve  the  transient  problem.   Baum  introduced  the  Singu- 
larity Expansion  Method  (SEM)  in  1971.   SEM  is  based  on  an 
extension  of  the  Heaviside  expansion  theorem,  used  in  circuit 
theory,  brought  into  the  realm  of  electromagnetic  fields. 
The  motivation  behind  most  work  in  SEM  has  been  the 
determination  of  EMP  threat  levels  in  military  electronic 
hardware.   SEM  represents  the  currents  induced  on  the 
scatterer  as  the  real  part  of  a  series  of  complex  exponen- 
tials such  that 

N 
f(t)  =  Re{   E   R   expCjs  t)}  (1) 

n=l   n       n 


where  R   is  the  amplitude  of  each  mode  of  a  complex  frequency 

s   [Ref.  1]  .   References  2  through  4-  address  current 
n  ° 


research  in  this  area. 
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In  many  cases,  such  as  in  radar  target  identification, 
the  surface  currents  on  the  target  are  of  importance  only  in 
the  context  of  generating  the  scattered  fields.   Another 
approach  to  the  problem  is  the  Finite  Element  Method  (FEM) 
described  in  Ref.  5.   This  method  offers  a  direct  solution 
of  the  scattered  electromagnetic  field  as  opposed  to  analy- 
sis of  surface  currents.   Although  the  approaches  are  dif- 
ferent, the  final  results  are  similar  since  the  scattered 
fields  are  related  to  the  surface  currents  by  Maxwell's 
equations . 

B.   GENERAL  GOALS 

In  order  to  study  transient  electromagnetics  and  to 
exercise  the  methods  previously  discussed,  facilities  must 
be  provided  for  experimentation.   With  this  idea,  the  Naval 
Postgraduate  School  Time  Domain  Measurements  Laboratory  was 
developed.   This  laboratory  provides  a  means  of  obtaining 
time  domain  measurements  of  the  transient  scattering  response 
of  an  object  to  incident  signals  produced  by  voltage  pulses 
or  steps  applied  to  a  transmitter  antenna  and  reflected  from 
the  target.   The  laboratory  was  constructed  primarily  to 
measure  the  transient  scattered  field  from  plane-symmetrical 
objects  for  impulse  response  determination  and  target 
classification.   This  is  consistent  with  the  long  range  goal 
which  is  to  provide  a  means  which  will  ultimately  establish 
the  feasibility  of  developing  advanced  radar  systems  which 
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are  capable  of  target  discrimination  and  classification  via 
time  domain  scattering  analysis.   In  addition,  research  in 
radar  cross  section  modeling  and  nuclear  EMP  vulnerability 
is  also  underway  which  extend  the  scope  of  the  original 
project . 

C.   FINAL  PHYSICAL  LABORATORY  DESCRIPTION 

Figure  1  shows  the  physical  configuration  of  the 
laboratory.   A  block  diagram  of  the  scattering  range  is 
shown  in  Figure  2.   The  signal  source  is  a  mercury  reed 
switch  type  pulse  generator  which  produces  a  250  picosecond 
risetime  voltage  pulse  of  calibrated  amplitude  and  variable 
pulse  width.   The  -incident  pulse  is  sent  through  50  ohm 
coaxial  cable  and  a  60  nanosecond  delay  line  to  a  twenty-one 
foot  monopole  transmitting  antenna  protruding  through  the 
aluminum  image  plane.   The  image  plane  isolates  the  trans- 
mitting and  receiving  antennas  from  the  electronic  equipment 
below.   The  incident  electromagnetic  field  radiates  outward 
from  the  transmitting  antenna.   This  field  illuminates  a 
target  on  the  image  plane,  and  subsequently,  the  target 
reradiates  a  scattered  electromagnetic  field.   The  incident 
field,  as  well  as  the  scattered  field,  are  received  by  the 
Transverse  Electric  and  Magnetic  (TEM)  horn  antenna.   The 
received  signal  is  routed  through  a  50  ohm  coaxial  cable  to 
a  digital  processing  oscilloscope  (DPO)  that  has  been 
triggered  by  a  portion  of  the  signal  from  the  pulse  genera- 
tor.  The  received  signal  is  sampled  in  the  DPO,  then 
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Figure  2.   Block  Diagram  of  Scattering  Laboratory 
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transferred  to  a  minicomputer  for  storage  and  further  signal 
processing.   Figure  3  shows  the  response  with  no  targets 
on  the  image  plane  to  a  50  volt  amplitude  pulse  of  20 
nanosecond  pulse  width  at  the  oscilloscope  before  any  signal 
processing.   This  compares  very  favorably  with  the  results 
in  reference  [6].   The  display  shows  a  20  nanosecond  time 
window  between  the  direct  transmission  and  the  erratic 
response  caused  by  the  end  of  the  incident  pulse  and 
reflections  from  the  tip  of  the  antenna.   This  time  window 
is  convenient  to  view  a  target  response  since  unwanted 
reflections  can  be  gated  out  in  time  with  the  oscilloscope. 
Figure  4  shows  the  response  to  the  same  input  pulse  with 
two  conically  shaped  metallic  targets  on  the  image  plane. 
The  target  response  is  seen  in  the  20  nanosecond  time 

window. 

* 

D.   MEASUREMENT  PROCEDURE 

The  response  in  Figure  4-  contains  the  superpositibned 
incident  electromagnetic  field  and  the  scattered  field  from 
a  target.   In  addition,  the  waveform  is  impregnated  with 
noise  from  internal  and  external  sources.   In  order  to 
determine  the  transient  response  of  a  target,  the  noise 
must  be  reduced  to  an  acceptable  level,  and  the  incident 
electromagnetic  field  must  be  subtracted  from  the  signal. 
In  any  sequence  of  measurements  of  the  transient  response, 
certain  procedures  must  be  carried  out.   First,  the  incident 
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Figure  3 . 
Response  to  50 
Volt  Amplitude. 
20  nanosec  Wide 
Pulse 


Figure  M- . 
Response  of  2 
Conical  Targets 
to  a  50  Volt, 
2  0  nanosec  Widef 
Incident  Pulse 
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signal  is  measured  without  a  target  on   the  image  plane. 
Results  similar  to  Figure  3  should  be  obtained.   N  successive 
incident  signals  are  transmitted  from  the  monopole  antenna, 
received  by  the  TEM  horn  antenna,  sampled  by  the  DPO ,  and 
transferred  to  the  minicomputer.   There  a  running  summation 
of  the  incoming  waveforms  is  tallied  until  all  of  the  N 
signals  have  been  processed.   A  time-average  of  the  sampled 
data  set  of  incident  waveforms  is  computed  by  the  equation 


1    N 
Al(k)  =  ±    E   A.(k)       k=l,2,...512  (2) 

N   .  ,   i 
1  =  1 


where  A. (k)  is  the  ith  waveform  transferred  from  the  digital 
i 

processing  oscilloscope.   The  value  of  N  is  made  large 
enough  to  sufficiently  reduce  the  noise  level,  as  will  be 
discussed  in  Chapter  II,  Section  C. 

Next,  a  target  is  placed  on  the  image  plane  at  a  parti- 
cular distance  and  aspect  angle  from  the  antennas.   Once 
again  N  successive  signals  are  processed  as  before,  then  a 
time-average,  Bl(k),  is  computed  using  equation  (2)  for  the 
sampled  data  set  of  waveforms  containing  both  the  incident 
signal  and  the  scattered  response.   These  waveforms  will  be 
referred  to  as  "augmented  waveforms". 

To  remove  the  incident  signal  from  the  augmented  wave- 
form, a  difference  C(k)  is  found  by 

C(k)  a  Bl(k)  -  Al(k)       k  =  l,2, ,512  (3) 
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C(k)  represents  the  sampled  target  response  convolved  with 
the  transient  responses  of  the  transmitting  and  receiving 
antennas . 
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II.   SCATTERING  LABORATORY  DEVELOPMENT 

A.   CONSTRUCTION  OF  PROTOTYPE  MODELS 

The  development  of  the  present  scattering  laboratory 
was  based  on  knowledge  attained  from  two  earlier  prototype 
laboratories  constructed  during  this  endeavor.   In  the 
development  of  the  first  model,  learning  the  principles  of 
scattering  range  operation  and  the  essentials  of  equipment 
interface  were  the  paramount  goals.   In  the  second  prototype, 
proper  antenna  selection  was  the  primary  consideration. 

The  first  model  consisted  of  a  22cm  x  22cm  aluminum 
image  plane  with  two  straight,  copper  wire  antennas  pro- 
truding through,  as  shown  in  Figure  5.   The  center  wire 
antenna  was  used  to  transmit  the  incident  signal  from  the 
pulse  generator.   The  receiving  antenna  was  placed  in  the 
corner  of  the  image  plane  and  connected  by  coaxial  cables 
to  the  sampler  head  of  the  sampling  oscilloscope.   Figure  7 
shows  the  response  to  a  37.5  volt,  incident  step  pulse  with 
a  10  nanosecond  pulse  width,  Figure  6.   In  reference  7,  a 
similar  configuration  is  described,  and  the  results  in  both 
cases  are  similar. 

The  initial  model  proved  that  the  electronic  equipment 
could  be  satisfactorily  interfaced  to  produce  an  inter- 
action between  the  two  antennas;  however,  several  factors 
inherent  with  this  model  made  quantitative  analysis  of  the 
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Figure  5 .   First  Prototype  Scattering  Range  Model 
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Figure  6 , 
37.5  Volt  10  nsec 
Pulse  Width, 
Incident  Pulse. 


Figure  7 . 
Response  to 
37. 5  Volt  Ampli- 
tude 10  nsec 
Pulse  Width 
Incident  Pulse 
(figure  6) 
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response  difficult.   These  factors  were  related  to  the  small 
size  of  the  image  plane  and  the  close  proximity  of  the 
antennas . 

To  overcome  these  problems,  a  larger  model  was  con- 
structed.  In  the  second  model,  a  longer  transmitting  antenna 
was  used,  and  two  broadband  receiving  antennas  were  tested 
to  improve  the  quality  of  the  response.   The  equipment  inter- 
face below  the  image  plane  was  essentially  identical  to  the 
first  model. 

The  length  of  the  transmitting  antenna  was  increased  to 
18  2  cm.   The  travel  time  for  a  pulse  to  reach  the  tip  of  the 
antenna  and  be  reflected  was  approximately  5  nanoseconds. 
This  provided  a  longer  time  window  on  the  oscilloscope  dis- 
play before  unwanted  reflections  appeared  in  the  signal. 

Two  receiving  antenna  designs  were  tested.   The  first 
antenna  was  a  coaxial  cone  antenna  (Figure  8).   Experimental 
and  theoretical  descriptions  are  contained  in  reference  8. 
The  received  response  of  this  antenna  to  an  incident  field 
is  a  differentiation  with  respect  to  time.   This  type  of 
response  is  useful  in  scattering  work  when  the  radiated  field 
is  a  step  function  in  time. 

Another  antenna  design  in  the  literature  is  the  TEM  horn 
antenna.   Analyses  in  the  frequency  and  time  domains  are 
found  in  references  9  and  10.   The  TEM  horn  antenna  can  be 
viewed  as  a  pyramidal  horn  with  open  sides.   See  Figure  9. 


22 


Figure  8 .   Coaxial  Cone  Antenna 
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Figure  9.   TEM  Horn  Receiving  Antenna 
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If  the  horn  flare  angle  9  and  the  plate  width  W  are  chosen 
properly,  a  constant  characteristic  impedance  can  be 
obtained. 

The  TEM  horn  antenna  was  selected  for  use  in  the 
scattering  laboratory.   Factors  influencing  the  choice  were 

1.  BETTER  DIRECTIVITY.   The  TEM  horn  antenna  can  be 
directed  at  a  target  anywhere  on  the  image  plane  to 
maximize  reception. 

2.  LESS  ANTENNUATION.   The  TEM  horn  antenna  had  less 
attenuation  than  the  coaxial  cone  antenna  tested. 

3.  MOBILITY  ON  THE  IMAGE  PLANE.   The  TEM  horn  antenna 
is  much  easier  to  move  around  on  the  image  plane 
due  to  its  design.   The  coaxial  cone  antenna  must 
remain  in  a  fixed  position  on  the  range  flush  with 
the  image  plane. 

4.  CHARACTERISTIC  IMPEDANCE.   The  characteristic  impe- 
dance of  the  TEM  horn  antenna  is  adjusted  by 
changing  the  flare  angle.   The  coaxial  cone  antenna 
requires  redesign  of  the  dimensions  to  alter  the 
characteristic  impedance. 

Having  made  the  receiving  antenna  selection,  Figure  10 
shows  the  final  configuration  of  the  second  model.   When  a 
3  5  volt  amplitude  pulse  with  a  20  nanosecond  pulse  width 
was  applied  to  the  transitting  antenna,  the  response  in 
Figure  11  was  obtained.   Although  the  results  did  not  match 
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Figure  11.   Second  Prototype  Scattering  Range  Model 
The  Response  to  a  35  Volt,  20  nanosec 
Pulse  Wide  Incident  Pulse. 
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those  in  Figure  3  and  in  reference  6,  there  were  fewer  oscil- 
lations than  in  the  first  model  indicating  the  TEM  horn  is 
a  satisfactory  receiving  antenna. 

B.   ANTENNA  DEVELOPMENT 
1.   Receiving  Antenna 

Having  made  the  preliminary  receiving  antenna 
selection  heuristically ,  a  theoretical  approach  was  pursued 
to  investigate  the  transmitting  and  receiving  characteristics 
of  the  TEM  horn  antenna  as  depicted  in  Figure  9 .   It  was 
assumed  that  the  flare  angle  0  and  the  plate  width  W  were 
chosen  so  that  the  characteristic  impedance  is  approximately 
constant  and  only  the  TEM  mode  propagates  in  the  horn.   Edge 
diffraction  and  fringing  fields  are  neglected.   The  field 
at  the  aperture  is  thus  given  by  the  TEM  wavefront 


V    -jkr 
E  (x,y,0)  =  -j2.  e     °  (4) 


where  V   is  the  amplitude  of  the  incident  voltage  at  the 
feedpoint,  h  is  the  height  of  the  aperture,  and  k  is  the 
wave  number  2tt/A.   The  E  field  closely  approximates  the 
spherical  wave  assumed  in  reference  C 10] .   Knowing  the 
field  at  the  aperture,   equivalent  magnetic  currents  are 
found  from  Harrington,  reference  11,  by 
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M  =  2E  x  n  (5) 

thus 

2V    -jkr 

M   =  — a  e     °  (6) 

x    h 

The  electric  vector  potential  F  is  given  by 


i  -jkr 

F  =  -±-  r     .  M ds  (7) 

4-tt  /  /    .        r 

aperture 


where  r  is  the  distance  from  the  aperture  to  an  observer 
on  the  image  plane  in  the  boresight  direction.   In  the  far- 
zone  of  the  aperture  F  then  becomes 

-jk(r+r  ) 
V  We       ° 

F   =  — 5= (8) 

x         2ttt 

The  far-zone  radiation  fields  form  a  TEM  wave  which  are 
related  to  the  electric  vector  potential  by 


H   =  -jwe   F  (9) 

X      J    O    X 


E   =  n  H   =  -jwe  n   F  (10) 

y    ox     J   o  o   x 


where  r\    =120tt  ohms  is  the  characteristic  impedance  of  free 


space.   The  radiated  fields  become 
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-jk(r+r    ) 

V    We                   ° 
tt  rad             o  /■    •         \  r  -,  i  \ 

H  = •    (-]we    )  (11) 

X  2TTT  O 


„     ,           -j60w£    V   We 
E  rad    = 2-5 (12) 


The  transmitting  transfer  function  is  given  in 
reference  9  by 


q  ffn  -  EradCf)    -J60  F(f)e"jkr  nq. 

V*'  '   V  (f)   "  [Z  (f)+Z  (f)Jr  {ld) 

o         o      g 


where  F(f)  is  the  radiation  function,  Z  (f)  is  the  input 

impedance  of  the  horn,  and  Zg(f)  is  the  generator  impedance 

Using  equation  12,  S"(f)  becomes 

-jk(r  +r) 
-  j  60  (27Tf)e  We     ° 
S  (f)  =  2 (11+) 


When  simplified,  (14)  becomes  in  terms  of  F(f) 

-  jkr 

F(f)  =   Z  (f)  +  Z  (f)   -27Tfs  We     °  (15) 

o      g  o 

Assuming  the  load  impedance  ZT  and  the  generator  impedance 
Z  "(g)  are  real  and  equal,  the  Rayleigh-Carson  reciprocity 
theorem  in  reference   12   provides  that 


F(f)  =  -h  (f)k  (16) 

e 
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By  using  (16)  the  effective  height,  h  (f),  becomes 

-jkr 
27Tf£  We     °(Z  (f)+Z  (f)) 

h  (f)  =  2 ^ S (17) 

e  k 


The  receiving  transfer  function,  S  (f),  can  be 
determined  with  the  aid  of  reference  9  as 


V  (f )      h  (f )Z. 

S  (f)  =  -Ji :-,  ve,f,  ^  (18) 

Elnc(f)     ZoCf)+ZL 


Substituting  equation  (17)  into  (18)  leads  to 

-jkr 

S  (f)  =  -ce  WzTe     °  (19) 

r         o   L 

Taking  the  inverse  Fourier  transform,  where 
s(t)  =  f"   S(f)ej2TTftdf  (20) 


the  receiving  transfer  function,  S  (f),  gives 


s  (t)  =  -ce  WzT5(t-r  /c)  (21) 

r         o   L     o 


while  the  transmitting  transfer  function,  S  (f ) ,  in  (14) 
yields  the  transmitting  impulse  response 


60s  W 

s  ft)  =- — =r2 —  -A.  5  (t_r  /c)  (22) 

t        r      dt       o 
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Comparing  (21)  and  (22)  shows  that  the  transmitting  impulse 
response  is  simply  the  time  derivative  of  the  receiving 
impulse  response.   This  is  in  agreement  with  the  results 
in  reference  [10] . 

Consider  equation  (21).   The  impulse  response  is  an 
impulse,  delayed  in  time  by  r  /c,  multiplied  by  a  negative 
constant  which  is  dependent  on  the  width  of  the  antenna 
aperture,  W,  and  the  load  impedance,  zT  .   This  theoretical 
result  can  be  accepted  if  it  can  also  be  shown  to  be  experi- 
mentally valid. 

The  far-zone  field  step  response  of  a  linear  dipole 
antenna  is  computed  in  references  7  and  13.   If  the  step 
response  of  a  dipole  antenna,  as  measured  by  the  voltage 
output  of  the  TEM  horn  antenna,  is  comparable  to  the 
theoretical  transmitted  field  results  in  references  7  and 
13,  then  two  implications  follow:   first,  the  TEM  horn  is 
apparently  in  the  far  field  of  the  dipole,  and. second,  the 
TEM  horn  produces  negligible  dispersion  and  distortion  of 
the  received  signal.   Thus  equation  (21)  can  be  verified. 

A  one  meter,  wire  monopole  antenna;  equivalent  to 
a  two  meter,  dipole  antenna;  was  constructed  on  the  image 
plane.   (The  7  meter  long  transmitting  monopole  antenna 
was  removed. )   The  meter  long  antenna  was  center-fed  with 
a  smoothed  step,  incident  signal  from  the  pulse  generator. 
Figure  12a  depicts  the  driving  voltage  measured  at  the 
antenna  feedpoint.   The  augmented  signal  in  Figure  12b  is 
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the  two  meter  dipole  response  measured  by  the  TEM  horn 
antenna  at  a  distance  of  270  cm  and  at  a  broadside  angle  of 
incidence.   Using  the  signal  processing  software  routine 
described  in  Chapter  II,  Section  E,  the  impulse  response  of 
the  dipole  antenna  was  computed  through  the  use  of  fast 
fourier  transforms  and  the  inverse  fourier  transform.   The 
impulse  response,  shown  in  Figure  13,  was  integrated  with  the 
use  of  a  signal  processing  read  only  memory  (ROM)  described 
in  Chapter  II,  Section  D.   Integration  of  the  impulse  res- 
ponse yields  the  step  response  of  the  coupled  antennas, 
Figure  1M-. 

This  experimentally  obtained  step  response  was  com- 
pared to  the  results  in  reference  13.   There  were  two 
factors  which  caused  the  experimental  and  theoretical 
results  to  differ.   First,  the  experimental  step  response 
lagged  the  theoretical  step  response  in  time;  second,   the 
experimental  voltage  and  a  reversed  polarity  compared  to 
the  theoretical  field.   Both  the  time  delay  and  the  negative 
factor  are  contained  in  equation  (21),  and  hence  are  attri- 
butable to  the  properties  of  the  TEM  horn  antenna.   Agreement 
between  the  two  results  is  quite  good.   See  Figure  15.   This 
validates  the  assumptions  used  to  establish  (21),  thus 
making  (21)  a  reasonable  representation  for  the  impulse 
response  of  the  TEM  horn  antenna. 
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2 .   Transmitting  Antenna 

The  implication  of  equation  (21)  is  that  a  receiving 
signal  waveform  will  be  guided  through  the  TEM  horn  antenna 
without  distortion.   Although  this  is  not  completely  true, 
since  there  will  be  some  losses,  dispersion,  and  distortion 
due  to  impedance  mismatches  and  multiple  modes,  this  assump- 
tion was  used  to  experimentally  determine  the  impulse  response 
of  the  transmitting  antenna.   By  disconnecting  the  coaxial 
cables  at  the  transmitting  and  receiving  antennas  and  re- 
connecting the  cables  as  shown  in  Figure  16 ,  the  input  signal 
x(t)  to  the  transmitting  antenna  can  be  measured  directly. 
Figure  17a  shows  the  incident  signal  to  be  a  smoothed  step 
function  with  approximately  a  250  picosecond  risetime.   The 
trailing  edge  of  the  pulse  was  time  windowed  out  using  the 
sampling  oscilloscope.   When  the  cables  were  reconnected  as 
in  Figure  1,  the  response  y(t)  was  measured  and  is  shown 
in  Figure  17b.   This  response  is  related  to  the  incident 
signal  by  the  double  time  convolution 


y(x)  =  [h(t)*s  (x)]*x(t) 
r 


where  h(t)  is  the  impulse  response  of  the  transmitting 

antenna,  s  (x)  is  the  receiving  impulse  response  of  the 

TEM  horn  antenna,  and  x=t-r  /c  is  the  time  delay.   Taking 

o 

the  magnitude  of  the  Fourier  transform  of  y(x)  yields 

|Y(f ) |  =  |H(f )S  (f) | |X(f ) |  (23) 
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Figure  16.   Coaxial  Cable  Connection  Used 
to  Measure  Incident  Pulse 
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The  magnitude  and  phase  plots  of  Y(f)  and  X(f)  are  shown  in 
Figure  18.   Dividing  these  results  gives 


G(f)|^|H(f)S  (f)|  =  |Y(f) |/|X(f) I  (24) 

i  i      r 


Taking  the  inverse  Fourier  transform  of  (24)  yields  an 
expression  g(t)  for  the  convolved  transient  impulse  responses 
of  the  transmitting  and  receiving  antennas  such  that 


g(t)  =  h(t)*s  (t)  (25) 

r 


Using  (21)  and  the  results  of  the  experimental  determination 
of  the  step  response  of  the  dipole,  g(t)  represents  the 
impulse  response  of  the  transmitting  antenna  such  that 


g(x)  =  -(ce  WzT )h(t)  (26) 

o   L 


where  T=t-r  /c.   Figure  19  shows  the  experimentally  deter- 
mined value  of  g(x). 

From  reference  14,  the  time  domain,  electromagnetic, 
far-zone  fields  transmitted  by  a  wire  antenna  are  given  by 


H(r,t)  =  -i—  /   |-  J   (r' ,T)xT^r  ds'  (27) 

4?rrc   s   8  t   s         r 


E(r,t)  =  ti  r  x  H(r,t)  "  (28) 
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where  r  is  the  distance  of  an  observer  from  the  antenna, 
Js(r',x)  is  the  surface  current  density  on  the  antenna  at 
a  distance  r'  from  the  origin,  and  t  is  the  propagation 
delay  time,  such  that  T=t-r/c.   From  equation  (27)  and  (28), 
the  far-zone  fields  radiated  are  proportional  to  the  time 
derivative  of  the  surface  current  density  on  the  antenna. 
Intuitively,  this  implies  that  the  transmitting  antenna 
functions  as  a  differentiator  with  respect  to  an  input 
incident  signal.   The  far-zone  field  response  to  a  step 
input  driving  voltage  will  be  impulse  shaped.   Similarly, 
if  an  impulse  driving  voltage  was  inputted,  the  far-zone 
field  response  would  be  the  derivative  of  an  impulse  or  a 
doublet  function.   This  is  the  result  that  would  be  anti- 
cipated for  the  impulse  response  of  the  transmitting  antenna, 
The  experimental  response  bears  this  supposition  out  in  that 
it  resembles  a  smoothed  doublet  function  as  shown  in  Figure 
19.   In  addition  to  this,  the  experimentally  determined 
impulse  response  can  be  shown  to  be  valid  by  comparison  with 
Figure  13 ,  which  is  the  impulse  response  of  the  one  meter 
wire  monopole  antenna  and  the  TEM  horn.   In  this  comparison, 
both  curves  are  identical  for  the  first  four  nanoseconds 
of  each  response.   After  this  time  period,  reflections  from 
the  end  of  the  one  meter,  monopole  antenna  make  the  impulse 
responses  differ.   In  Figure  19,  the  transmitting  antenna 
was  long  enough  so  that  reflections  from  the  end  of  the 
antenna  were  time  windowed  out  by  the  oscilloscope;  hence, 
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the  impulse  response  contains  no  reflections.   Since  the 
step  response  of  the  one  meter  wire  antenna  was  proved  to 
be  valid,  the  impulse  response  of  the  one  meter  antenna 
is  likewise  correct.   From  this  inference,  a  comparison 
of  the  similarities  in  Figures  13  and  19  make  the  impulse 
response  of  the  transmitting  antenna  a  valid  representation 
for  the  first  four  nanoseconds  of  time. 

There  are  several  constraints  imposed  by  the  TEM  horn 
antenna  which  may  cause  Figure  19  to  differ  from  its  actual 
value.   The  first  limitation  is  due  to  the  TEM  horn  antenna 
not  being  in  the  far-zone  field  of  the  transmitting  antenna; 
and  the  second,  which  is  a  consequence  of  the  first,  is  than 
the  flare  angle  of  the  TEM  horn  is  not  wide  enough  for  the 
horn  to  receive  the  electromagnetic  fields  at  its  aperture 
when  the  fields  are  radiated  by  the  upper  portions  of  the 
transmitting  antenna.   A  third  constraint  is  based  on  the 
assumption  that  a  TEM  wavefront  exists  at  the  aperture.   This 
is  valid  only  as  long  as  the  wavelengths  of  the  received 
signal  are  large  compared  to  the  TEM  horn  aperture  dimen- 
sions.  There  is  a  possibility  that  higher  frequencies  may 
be  distorted  or  attenuated  while  propagating  through  the 
horn  antenna.   This  would  account  for  smoothing  which  may 
occur  in  the  measured  scattered  response  from  a  target  on 
the  image  plane. 
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C.   SIGNAL  PROCESSING  SCHEME 

A  significant  problem  in  signal  processing  schemes  is 
the  estimation  of  signal  parameters  such  as  the  mean,  root 
mean  square,  and  variance  in  the  presence  of  noise.   The 
noise,  which  appears  in  the  measurements  made  in  the 
scattering  laboratory,  is  generated  by  external  and  inter- 
nal sources.   The  broadband  receiving  characteristics  of 
the  TEM  horn  antenna  make  the  system  especially  susceptible 
to  both  man-made  and  natural  external  noise  sources.   Man- 
made  sources  include  television,  radio,  and  radar  emissions 
from  the  local  area.   Natural  external  sources,  dependent  on 
atmospheric  conditions,  change  throughout  the  day.   Thermal 
agitation  of  electrons  is  increased  during  mornings  and 
afternoons  when  the  antennas  and  the  image  plane  are 
directly  exposed  to  the  sun.   Conversely,  thermal  excita- 
tion of  electrons  is  reduced  at  night  when  moisture  con- 
denses on  the  image  plane  and  cooling  of  surface  temperatures 
results . 

Internal  noise  has  several  sources  within  the  system. 
Mutual  interactions  between  the  transmitting  and  receiving 
antennas,  as  well  as  the  interaction  between  the  antennas 
and  a  target,  produce  invalid  signal  responses  unless 
sufficient  separation  is  maintained.   To  minimize  antenna 
interactions,  the  TEM  horn  antenna  was  placed  approximately 
three  meters  from  the  transmitting  antenna.   With  power 
reduction  inversely  proportional  to  the  square  of  the 
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distance  between  the  antennas ,  the  interactions  are  reduced 
by  a  factor  of  about  10  as  compared  to  a  separation  of  one 
meter  or  less  between  the  antennas  used  in  the  prototype 
models.   Separation  between  the  antennas  and  a  target  is 
dependent  on  the  aspect  angle  of  the  target  to  be  viewed. 
When  placing  the  target  on  the  image  plane,  it  should  be 
positioned  as  far  away  from  the  antennas  as  possible  while 
still  receiving  a  target  response  of  sufficient  amplitude. 
In  this  way  the  measurement  of  far-field  scattering 
quantities  Is  optimally  approximated. 

The  antenna  responses  also  include  pollution  by  internal 
noise  sources.   Multimode  propagation,  edge  diffraction,  and 
fringing  fields  cause  distortion  and  dispersion  of  the 
received  signal  which  alter  the  TEM  wavefront  at  the  apera- 
ture  of  the  horn  antenna. 

Another  significant  internal  noise  source  in  the  system 
is  the  pulse  generator.   Johnson  noise,  caused  by  circuit 
components  in  the  pulse  generator,  is  always  present. 
Another  problem  is  created  when  the  pulse  generator  is 
interfaced  with  the  digital  processing  oscilloscope.   When 
the  computer  program  subtracts  the  incident  averaged 
waveform  from  the  augmented  averaged  waveform,  there  is 
sometimes  a  shift  between  the  two  signals  on  the  horizontal 
time  axis.   This  shift  produces  a  noise  spike  in  the 
resulting  target  response.   The  amplitude  of  this  spike  is 
random,  but  the  spike  occurs  at  the  same  points  in  time 
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as  the  leading  and  falling  edges  of  the  incident  signal. 
These  noise  spikes  introduce  erroneous  high  frequency 
components  in  the  spectrum  and  thus  invalidate  the  results 
To  illustrate  this  phenomenon,  a  set  of  fifty  signals  was 
measured  and  averaged.   Without  changing  any  settings  on 
the  oscilloscope  or  the  pulse  generator,  another  set  of 
fifty  signals  was  measured  and  averaged.   The  two  averaged 
signals  are  shown  in  Figure  20.   The  difference  between 
the  two  waveforms  is  shown  in  Figure  21.   The  amplitudes 
of  the  noise  spikes  are  significantly  greater  than  the 
other  noise  present  in  the  system.   Signal  averaging  does 
not  decrease  the  noise  spike.   Increasing  the  number  of 
signals  to  be  averaged  increases  the  time  to  measure  a  set 
of  waveforms.   This  often  causes  the  shift  on  the  time 
axis  to  drift  further,  resulting  in  an  even  larger  noise 
spike.   Since  the  amplitude  of  the  spike  is  random,  better 
results  can  be  usually  obtained  by  repeating  the  complete 
set  of  measurements. 

In  addition  to  the  waveforms  displayed  in  Figures  2  0 
and  21,  the  peak  values,  root  mean  square  and  mean  values 
are  also  given.   Even  though  the  incident  and  augmented 
waveforms  in  Figure  2  0  are  two  sets  of  waveforms  from  the 
same  ensemble,  the  values  of  the  statistical  parameters 
differ  to  some  degree.   In  order  to  make  the  estimate  of 
the  mean  an  accurate  approximation  of  the  true  mean, 
ensemble  waveform  averaging  is  used  to  reduce  noise  in  the 
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system.   To  determine  the  effectiveness  of  waveform  averaging, 
the  following  tests  were  made.   First,  the  pulse  generator 
was  disconnected  from  the  coaxial  cable  leading  to  the 
transmitting  antenna.   In  place  of  the  pulse  generator, 
a  matched  5  0  ohm  load  termination  was  connected  to  the 
cable.   The  image  plane  was  cleared  of  targets,  and  the  pulse 
generator  was  used  only  to  trigger  the  oscilloscope.   By 
measuring  the  received  signals  with  the  TEM  horn  antenna, 
the  amount  of  noise  in  the  system  could  be  determined.   If 
an  incident  noise  waveform  and  augmented  noise  waveform 
were  both  measured  and  then  subtracted,  the  resulting  noise 
is  what  could  be  expected  to  be  present  in  the  target  res- 
ponse.  Three  different  values  of  N,  the  number  of  waveforms 
to  be  averaged,  were  selected;  these  were  N=l,  N=10,  and 
N=50.   Figures  22  and  23  show  the  noise  present  for  N=l  and 
N=5  0  cases.   The  following  chart  shows  the  effect  that 
increasing  N  has  on  noise  reduction  for  the  target  response. 
These  results  compare  very  favorably  with  the  theoretical 
reduction  of  the  root  mean  square  value  by  the  square  root 
of  N  produced  by  ensemble  averaging  of  random  variables.   It 
was  determined  that  N=50  reduced  the  noise  to  an  acceptable 
level . 
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N=l       N=10  N=50 

Maximum  value  (volts)     0.0131     0.0048  0.0017 

Minimum  value  (volts)     -0.0141    -0.0041  -0.0019 

Rms  value  (volts)         0.0035     0.0011  0.0005 

Mean  value  (volts)        0.0002     0.0001  0.0000 

In  the  computation  of  the  statistical  parameters,  the 

A 

estimate  of  the  mean  value  X  is  computed  by 


512 
X  =  (1/512)    I  X(i)  (29) 

i=l 


where  X(i)  is  the  ith  sample  of  the  512  samples  made  by  the" 
digital  processing  oscilloscope.   Use  of  equation  (29)  assumes 
ergodicity,  that  is  a  time  average  is  equivalent  to  the 
ensemble  average  of  the  waveforms. 

The  estimate  of  the  root  mean  square  value  X(RMS)  is 
given  by 


/  512         ? 

X(RMS)  =/  (1/512)   Z    CX(i)]  (30) 

i  =  l 


An  estimate  of  the  variance  is  not  computed  in  the  program 
but  is  easily  determined  from  the  statistical  parameters 
already  computed.   The  estimate  of  the  variance  Var(X) 
is  determined  by 
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512  2 

Var(X)  =  (1/511)   E   Cx(i)  -  Xl  (31) 

i  =  l 


where  the  value  of  511  or  N-l  is  an  unbiased  estimator. 

D.   DATA  PROCESSING  HARDWARE 

1 .   Digital  Processing  Oscilloscope 

The  digital  processing  oscilloscope  (DPO)  used  in 
the  scattering  range  is  arranged  to  make  sampling  measure- 
ments of  fast  waveforms.   In  addition  to  this,  the  DPO  is 
capable  of  storing  four  waveforms  in  an  integral  memory  and 
is  also  a  computer  compatible  oscilloscope  which  provides 
an  interface  to  the  system  minicomputer. 

The  digital  processor  unit  of  the  DPO  consists  of 
six  sections  and  an  asynchronous  bus.   The  system  bus 
operates  on  a  first-come,  first-served  basis  and  provides 
all  sections  of  the  processor  with  parallel  access  to 
power,  address  lines,  control  lines,  and  data  lines.   Figure 
24  shows  a  block  diagram  of  the  processor  unit. 

Incoming  signals  to  the  DPO  are  received  by  the 
acquisition  unit  of  the  oscilloscope.   Signal  flow  between 
the  CRT  display,  the  acquisition  unit,  and  the  processor  is 
then  handled  by  the  signal  interface  unit.   Signal  sampling 
for  the  processor  is  done  by  the  sample  and  hold  section. 
Horizontal  and  vertical  samples  are  made  every  5.5  micro- 
seconds regardless  of  oscilloscope  sweep  speed.   For  sweeps 
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slower  than  0.5  milliseconds/division,  all  512  horizontal 
sample  points  are  filled  in  a  single  sweep.   When  repetitive 
waveforms  are  sampled  at  sweep  speeds  faster  than  0.5 
milliseconds/division,  updated  samples  are  continuously 
provided  while  the  DPO  is  in  the  STORE  mode.   When  the 
processor  is  placed  in  the  HOLD  mode,  a  complete  waveform 
is  held  in  memory  and  is  composed   of  old  and  new  samples. 
For  non-repetitive  waveforms  at  sweep  speeds  faster  than 
0.5  milliseconds/division  not  all  512  horizontal  memory 
locations  are  filled,  and  resolution  of  the  waveform  is 
reduced.   Graphs  are  provided  in  the  users  manual  which 
give  the  time  required  to  fill  the  memory  for  various 
sweep  speeds. 

Conversion  of  the  analog  signal  to  a  digital  signal 
compatible  with  memory  storage  is  accomplished  by  the 
analog  to  digital  (A-D)  converter.   The  horizontal  sample 
is  converted  to  a  nine  bit  word  (2**9  =512  points),  and 
the  vertical  sample  is  converted  to  a  ten  bit  word  (2**10  = 
1024  points) . 

After  the  analog  signals  are  digitized,  the  memory 
section  stores  the  waveforms  as  512  horizontal  points  with 
1024  possible  levels  of  vertical  resolution  at  each  point. 
Each  of  the  512  horizontal  locations  contains  sequential 
samples  which  are  single-valued.   The  memory  consists  of 
4096  words  partitioned  into  eight  distinct  divisions  of 
512  words  each.   Four  divisions  are  reserved  for  the  storage 
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of  waveforms.   Each  of  the  four  divisions  has  four  fields 
associated  with  it.   These  fields  are  used  for  the  storage 
of  waveform  information  such  as  scale  factor  data  or  mini- 
computer generated  messages.   Information  stored  in  the  four 
fields  can  be  displayed  simultaneously  with  the  waveform 
on  the  CRT  display. 

The  display  generator  of  the  processor  unit  is  used 
to  display  waveform  data  on  the  CRT  which  has  been  placed 
in  memory  by  the  A-D  converter  or  the  minicomputer.   Two 
modes  of  CRT  display  are  possible,  an  X-Y  mode  or  a  Y-Time 
mode.   Only  the  Y-Time  mode  is  compatible  with  memory  storage 
When  a  sampled  waveform  is  displayed,  the  display  generator 
uses  a  linear  interpolation  feature  to  connect  all  adjacent 
points  with  line  segments. 

Finally,  the  Input/Output  interface  section  of  the 
processor  acts  as  a  bilateral  link  between  the  processor 
and  minicomputer;  the  front  panel  section  allows  operation 
of  the  DPO  from  the  front  panel ;  and  the  readout  interface 
converts  readout  information  from  the  acquisition  unit  into 
ASCII  code  for  storage  in  memory  and  reconverts  it  into 
readout  information  for  the  display  unit. 
2 .   Minicomputer 

Control  of  data  acquisition  and  signal  processing 
is  managed  by  a  64k  minicomputer  system  controller.   Figure 
2  5  presents  a  block  diagram  of  the  controller  hardware 
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which  is  interfaced  with  the  DPO  and  a  hard  copy  display- 
unit.   The  keyboard  is  the  primary  input  device  for  the 
system.   Each  time  a  key  is  pressed  an  ASCII  coded  represen- 
tation is  written  into  a  random  access  memory  (RAM)  location 
called  a  line  buffer  and  displayed  on  the  CRT.   When  the 
return  key  is  pressed,  the  contents  of  the  line  buffer  is 
written  into  memory  by  the  minicomputer  processor  provided 
the  input  conforms  to  the  proper  syntax.   If  not,  the  message 
is  displayed  again  on  the  CRT  for  editing  along  with  a  syntax 
error  cursor  indicating  the  location  of  the  error. 

The  two  magnetic  tape  units  each  provide  300k  bytes 
of  additional  data  storage  capability  for  the  system.   The 
information  storage  is  permanent  and  may  be  a  software  pro- 
gram or  data  to  be  processed  by  a  program.   All  data  flow 
between  the  magnetic  tape  unit  and  the  RAM  must  pass  through 
the  processor. 

The  RAM  provides  5  6k  of  temporary  storage  memory  for 
keyboard  entries,  program  instruction,  and  intermediate  pro- 
cessing results.   The  contents  of  the  RAM  are  destroyed 
when  written  over  or  power  is  removed  from  the  system. 

The  system  processor  is  the  routing  and  computing 
device  for  the  minicomputer.   The  processor  directs  all 
system  operations,  decodes  program  instructions,  and  performs 
arithmetic  and  logic  operations.   An  integral  part  of  the 
processor  is  the  firmware  permanently  fixed  in  the  system 
which  interprets  program  instructions  for  processing.   This 
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is  known  as  the  BASIC  interpreter.   Also,  added  are  three 
signal  processing  read  only  memory  (ROM)  packs  of  up  to 
16k  bytes  of  specialized  programs.   These  ROM  packs  perform 
such  operations  as  fast  Fourier  transforms  (FFT),  inverse 
fourier  transforms  (IFT),  convolution,  correlation,  and  high 
level  text  editing  operations.   The  programming  language  used 
by  the  processor  is  Beginner's  Ail-Purpose  Symbolic  Instruc- 
tion Code  (BASIC)  developed  in  the  nineteen  sixties  at 
Dartmouth  College.   Extensions  of  this  language  are  incor- 
porated for  graphics  displays  and  for  transfer  of  data 
between  internal  and  external  devices. 

The  graphic  system  display  is  a  cathode  ray  tube 
storage  device  which  allows  for  interaction  between  the 
processor  and  the  operator  through  the  keyboard.   The  dis- 
play screen  retains  the  written  image  for  viewing  and  hard 
copy  reproduction.   The  addressable  points  on  the  CRT  are 
defined  in  terms  of  graphic  display  units  (GDU).   There 
are  130  horizontal  addressable  GDU ' s  and  100  vertical 
addressable  GDU's. 

The  general  purpose  interface  bus  (GPIB)  links  all 
peripheral  devices  to  the  system.   GPIB  specifications  con- 
form to  IEEE  standard  M-88-1975.   The  transfer  mode  between 
devices  is  byte-serial,  bit-parallel.   External  devices 
are  serviced  by  interrupt  instructions.   Enable,  polling 
and  data  transfer  commands  are  available  under  program 
control . 
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E.   SOFTWARE  DEVELOPMENT 

The  following  logic  flow  diagram  presents  an  overview 
of  the  software  program  developed  to  perform  all  signal 
processing  necessary  from  the  initial  acquisition  of  a 
signal  from  the  receiving  antenna  on  the  image  plane  to  a 
graphical  display  of  the  impulse  response  of  a  target. 
Appendix  A  contains  the  software  program 
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Figure  26.   Logic  Flow  Diagram 
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The  program  is  designed  to  request  N,  the  number  of 
waveforms  to  be  averaged,  interactively  from  the  operator 
through  the  minicomputer  keyboard.   After  the  value  of  N 
has  been  entered,  the  program  returns  that  value  for  veri- 
fication and  asks  the  operator  if  the  image  plane  is  pro- 
perly set  up.   After  receiving  an  affirmative  response 
of  "YES"  from  the  operator,  the  program  determines  if  the 
signals  to  be  averaged  will  contain  a  target  response 
(augmented  waveform)  or  if  the  range  will  be  clear  of  a 
target  and  only  the  direct  response  (incident  waveform) 
will  be  measured.   This  requires  the  operator  to  enter  a 
value  K  of  "1"  for  incident  waveform  averaging  or  "2"  for 
augmented  waveform  averaging.   When  this  value  K  is 
entered,  the  program  returns  that  value  for  verification. 
The  operator  then  inputs  the  command  "GO"  and  operation 
of  the  DPO  is  taken  over  by  the  minicomputer. 

Lines  280-410  are  a  FOR/NEXT  loop  which  instruct  the 
DPO  to  store  then  hold  a  waveform  in  memory.   Lines  310- 
315  from  a  BUSY/WAIT  loop  which  slows  the  execution  speed 
of  the  program  to  allow  for  waveform  storage  time.   The 
minicomputer  processor  then  prepares  to  transfer  the  digi- 
tized waveform  from  the  DPO  memory  to  the  minicomputer 
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memory.   After  the  interface  bus  is  clear,  the  transfer 
is  made  on  command.   Once  the  waveform  is  in  the  mini- 
computer memory,  the  process  is  repeated  until  N  waveforms 
have  been  acquired.   Each  time  a  waveform  is  transferred 
from  the  DPQ  ,  it  is  added  to  the  sum  of  the  previous 
waveforms  until  the  resultant  matrix  of  time-sampled  values 
B  is  obtained  where 

N 
B  =  Z   A.  (32) 

i=l   X 

For  incident  waveforms  the  ensemble  averaged  time-sampling 
matrix  Al  is  determined  where 

Al  =  B/N  (33) 

For  augmented  waveforms  the  ensemble  averaged  time-sampling 
matrix  Bl  is  computed  such  that 

Bl  =  B/N  (34) 

The  program  determines  next  if  all  signal  averaging  is 
complete.   If  not,  the  value  of  K  is  changed  and  the  program 
branches  back  to  the  FOR/NEXT  loop  to  process  another  set  of 
waveforms.   When  the  averaging  is  complete,  the  incident 
waveform  is  subtracted  from  the  augmented  waveform  such  that 

C  =  Bl-Al  (35) 
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Lines  700-940  calculate  the  maximum  and  minimum  values, 
mean  values,  and  root  mean  square  values  of  waveforms  C,  Bl, 
and  Al .   The  maximum  and  minimum  values  are  used  for 
plotting  purposes  in  the  program  to  determine  the  vertical 
window  length  of  the  graphical  displays. 

The  first  graphical  results  present  the  incident  and 
augmented  averaged  waveforms  along  with  the  peak  values, 
mean  and  rms  values  of  Bl  and  Al.   Refer  to  Figure  2  7  for 
an  example  of  this  graphical  format.   Two  subroutines  are 
used  to  scale  the  sampled  waveform  data  and  to  graphically 
plot  the  waveforms . 

Before  subroutine  4000  is  called,  each  of  the  512  hori- 
zontal points  has  some  discrete  value  between  0  and  1024. 
Subroutine  4000  transfers  the  voltage  scale  factor  and  the 
time  scale  factor  from  field  (F0)  commencing  at  address  2  048 
of  waveform  A  in  the  DPO  to  the  minicomputer  processor.   The 
sampled  data  values  are  then  multipled  by  the  appropriate 
scale  factors . 

Subroutine  5000  determines  the  scale  factor  of  the 
ordinate  of  each  graph  based  on  the  maximum  and  minimum 
values  previously  determined  then  plots  a  smooth  curve  for 
each  waveform.   In  a  similar  manner  the  target  response  C 
is  scaled,  then  plotted  using  subroutine  5000.   This  graph 
is  followed  by  the  maximum,  minimum,  mean  and  rms  values 
of  the  target  response.   Figure  2  8  provides  an  example. 
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The  target  response  C  represents  the  convolved  tansient 
responses  of  the  transmitting  and  receiving  antennas  con- 
volved along  with  the  transient  response  of  the  target.   In 
order  to  isolate  the  transient  target  response,  additional 
signal  processing  is  necessary.   First,  the  fast  Fourier 
transform  (FFT)  of  waveforms  C  and  Al  are  computed  with 
the  use  of  one  of  the  signal  processing  ROM's.   The  FFT 
command  performs  a  fast  calculation  of  the  discrete  Fourier 
transform  (DFT)  given  by 


C(k)    =         E         c(n)e"j27Tnk/572      k=0,l,2, 255  (36) 

n  =  0 


where  C(k)  is  the  kth  Fourier  coefficient  and  c(n)  refers 
to  the  (n+l)th  element  of  the  real  sampled  data  of  waveform 
C.   A  similar  set  of  values  is  generated  for  waveform  Al . 
In  this  particular  FFT  algorithm,  C(l)  and  AMD  represent 
the  direct  current  value  of  those  signals.   C(2)  and  AM  2) 
represent  the  value  of  the  FFT  at  the  Nyquist  frequency  which 
is  one-half  the  sampling  frequency.   The  remaining  values 
C(3)-C(51D  and  A1(3)-AM51D  represent  the  interleaved  real 
and  imaginary  terms  for  the  Fourier  coefficients. 

The  "POLAR"  command  allows  for  easier  interpretation  of 
the  FFT  results  by  converting  the  rectangular  Fourier 
coefficients  into  polar  form  where  array  components  are 
expressed  in  terms  of  magnitude  and  phase  angle.   Lines 
1930-22M-0  generate  graphical  displays  of  magnitudes  and 
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phases  of  the  FFT's  of  waveforms  C  and  Al .   See  Figure  30. 
Since  the  time  domain  waveforms  are  real,  the  FFT's  will  be 
even  and  symmetric  about  the  vertical  axis;  thus  only  256 
points  are  required  for  the  magnitude  of  the  FFT  and  256 
points  for  the  phase  of  the  FFT. 

At  this  point  the  FFT  of  waveform  Al  is  the  product  of 
the  system  transfer  functions  of  the  transmitting  and 
receiving  antennas  such  that 

Al(k)  =  H(k)  Sr(k)      k=0,l, 255  (37) 

where  H(k)  is  the  kth  component  of  the  transfer  function  of 
the  transmitting  antenna  and  Sr(k)  is  the  kth  component 
of  the  transfer  function  of  the  receiving  antenna.   The  FFT 
of  waveform  C  is  the  product  of  the  transfer  functions  of 
the  transmitting  antenna,  receiving  antenna,  and  the  target 
thus 

C(k)  =  T(k)  H(k)  Sr(k)  k=0,l, 255        (38) 

where  T(k)  is  the  kth  component  of  the  target  transfer  func- 
tion.  The  isolate  T(k),  the  magnitude  of  (38)  is  divided 
by  the  magnitude  of  (37)  so  that 

|T(k)|   =   |C(k) |/|Al(k) |       k=0,l, 255        (39) 
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and  the  phase  of  Al(k)  is  subtracted  from  the  phase  of  C(k) 
such  that 

<T(k)  =  <C(k)  -  <Al(k)        k=0,  255         (40) 

Having  determined  the  magnitude  and  phase  of  T(k),  it  is  a 
simple  matter  ot  convert  TCk)  from  polar  form  to  rectangular 
form.   This  is  done  in  lines  2340-2370  of  the  program.   Fin- 
ally, the  inverse  Fourier  transform  (IFT)  of  TCk)  is  calcu- 
lated by 

512 
h  Cn)  =  (1/512)    E    TCk)  e  +  32Tmk/512   n=0,l,..511 
t  k=0 

(41) 

where  h  (n)  is  the  (n+l)th  element  of  the  sampled  data 

impulse  response  of  the  target.   h  (n)  is  plotted  as  a 

function  of  time  by  lines  2430-2540  in  the  program. 

To  evaluate  the  signal  processing  scheme  of  the  program, 

an  ideal,  time  domain  step  function  described  by 


1.0     k=0,l, 100 

X(k)  =  {0       k=101,102. . .511  (42) 


was  input  into  the  program  as  the  averaged  incident,  sampled 
data  waveform;  and  the  averaged  augmented,  sampled  data 
waveform  was  input  as 

k       k=0,l, 100 

y(k)  =  {  200-k    k=101,102, 200 

0        k-201,202, 511  (43) 
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See  Figure  28.   When  the  program  was  executed,  the  target 
response  was  given  as  a  triangular  waveform  shown  in  Figure 
29.   This  is  analogous  to  imputting  a  square  wave  into  an 
unknown  linear  system  obtaining  the  response  as  shown  in 
Figure  27. 


tL 


x(k) 


-S 


h(k) 


A. 


-> 


y(k) 


Figure  2  7 

y(k)  is  expressed  mathematically  as  the  convolution  of 
x(k)  and  h(k)  such  that 

511 
y(k)  =  x(k)-h(k)  =   Z    h(i)  x(k-i)   for  k-0,l,...511 

i  =  0 

(44) 

Since  x(k)  and  y(k)  are  known,  h(k),  the  system  impulse 

response  needs  to  be  determined.   This  is  precisely  the 

main  objective  of  the  program. 

Given  y(k)  and  x(k),  the  FFT's  of  these  functions  are 

computed  and  displayed  by  the  program  as  shown  in  Figure  3  0 

X(nf ) ,  the  nth  component  of  the  Fourier  transform  of  x(k), 

has   a    [sin(2imf )  ]/(2irnf )    spectrum,    and   Y(nf )  ,    the   nth 
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component  of  the  Fourier  transform  of  y(k),  has  a 
[sin(2  nf )/(2-rrnf )]**2  spectrum.   H(nf)  is  then  found  by  (39) 
and  C+0),  and  the  resulting  spectrum  has  a  [sin(2irnf )  ]/ (2irnf ) 
shape . 

After  the  IFT  of  H(nf)  is  computed,  the  impulse  response 
is  obtained  and  shown  in  Figure  31.   This  demonstrates  the 
validity  of  the  program  in  finding  the  impulse  response  of  a 
system.   By  applying  this  same  technique  to  the  waveform 
signals  obtained  on  the  scattering  range,  the  impulse  response 
of  a  target  on  the  image  plane  can  be  determined. 

F.   TARGET  FABRICATION 

The  targets  fabricated  for  use  on  the  image  plane  con- 
sist of  several  axisymmetric  configurations ,  where  some  are 
made  of  aluminum  while  others  are  made  of  plexiglas  dielec- 
tric material.   Figure  32  shows  a  sample  of  the  target 
inventory.   When  placed  on  the  scattering  range,  the  targets 
are  "mirror-imaged"  by  the  surface  of  the  image  plane.   In 
this  manner  a  hemispherical  target  produces  the  scattered 
response  of  a  sphere  due  to  reflections  from  the  conducting 
surface  of  the  image  plane . 

Three  basic  designs,  a  cone,  a  cylinder,  and  a  hemisphere, 
were  machined  from  three  inch  diameter  and  six  inch  diameter 
aluminum  stock.   Later,  some  of  the  same  designs  were 
machined  from  a  plexiglas  material.   Appendix  B  contains 
the  specifications  of  the  dielectric  material.   The  dimensions 


Figure  32.   Sample  of  Representative  Targets 
Used  on  the  Image  Plane 
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of  the  plexiglas  targets  were  matched  to  those  of  the  alumi- 
num targets  so  comparisons  between  the  conducting  properties 
of  the  two  materials  could  be  made.   In  addition  to  these 
targets,  a  bisected  version  of  the  cylinder  was  fabricated  to 
increase  the  number  of  aspect  angles  that  this  target  could 
be  viewed  and  to  change  the  polarization  of  the  incident 
fields.   See  the  foreground  of  Figure  32.   Several  targets 
that  were  easy  to  fabricate  have  also  been  added  to  the 
inventory.   These  consist  of  two  aluminum,  rectangular  plates 
and  several  metallic  wires  of  different  lengths. 

Target  resolution  in  scattering  measurements  is  a  problem 
dependent  on  the  dimensions  of  the  target.   Several  of  the 
smaller  targets,  such  as  the  three  inch  diameter  hemisphere, 
have  scattering  responses  which  produce  signal  to  noise 
ratios  low  enough  to  be  indistinguishable  from  noise  in  the 
system  regardless  of  their  position  on  the  image  plane. 
Even  signal  averaging,  used  to  reduce  noise  levels,  is  in- 
effective in  enhancing  the  response. 

In  general,  there  is  no  rule  of  thumb  on  which  to  base 
the  criteria  for  dimension  designs.   This  is  dependent  on  the 
type  of  scattering  response  produced  by  a  target  and  on  the 
noise  present  in  the  system.   As  an  example,  a  wire  has  a 
damped  sinusoidal  target  response,  and  usually  can  be 
detected  by  the  DPO  even  with  a  length  dimension  slightly 
less  than  twelve  inches.   A  target  like  the  sphere,  whose 
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time  domain  response  is  fast,  is  very  difficult  to  detect. 
The  three  inch  diameter  sphere  has  never  been  detected  on 
the  image  plane  regardless  of  target  position. 

Solutions  to  this  problem  are  addressed  in  Chapter  IV, 
CONCLUSIONS  AND  RECOMMENDATIONS.   In  addition  to  the  obvious 
solution  of  building  larger  targets,  use  of  a  higher  voltage 
pulse  generator,  as  planned,  may  improve  target  resolution. 
With  a  1000  volt  signal  amplitude,  100  picosecond  risetime, 
the  scattered  responses  of  the  smaller  targets  may  be  in- 
creased enough  to  be  identifiable.   This  remains  to  be 
tested. 
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III.   SCATTERING  LABORATORY  MEASUREMENT  RESULTS 

Throughout  much  of  the  literature  on  time  domain  scat- 
tering as  in  references  14  and  15  for  example,  a  regulari- 
zation  of  an  impulse  is  used  as  the  incident  waveform  to 
illuminate  targets.   This  incident  field  is  referred  to  as 
a  smoothed  impulse  and  has  the  Gaussian  form 

2         2 
-a   t-t 

E(t)  =  E  e       max  (45) 

o 


where  E   is  the  peak  amplitude  of  the  incident  field, 

t     is  the  time  of  peak  amplitude  or  the  signal  mean,  and 
max 

a  is  a  spread  parameter  related  to  the  standard  deviation 
a   of  the  Gaussian  curve  by 


2     1 
az  =  -±j  (i+6) 

2a 


This  smoothed  impulse  is  a  convenient  waveform  because  it 
has  rapidly  decaying  time  and  frequency  domain  represen- 
tations, and  it  is  well-suited  for  numerical  solution  of  the 
time  domain  scattering  interaction. 

If  the  incident  field  radiated  from  the  transmitting 
antenna  is  approximated  by  this  Gaussian  waveform,  the 
smoothed  impulse  response  of  a  target  can  be  measured  by 
subtracting  the  incident  waveform  from  the  augmented  waveform, 
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In  this  manner  the  "target  response"  determined  by  the 
software  routine  is  the  smoothed  impulse  response  of  a 
target.   This  result  can  then  be  compared  to  theoretical 
results  in  order  to  ascertain  the  validity  of  scattering 
range  measurements. 

A  thin  wire  scatterer  with  a  50  cm  length  and  a  0.161 
cm  diameter  was  selected  as  the  target.   The  scattering  pro- 
perties of  a  wire  are  understood,  and  the  numerical  solution 
of  its  scattered  fields  is  well-documented. 

The  experimental  measurements  were  obtained  as  follows. 
Without  the  target  on  the  image  plane,  5  0  incident  wave- 
forms were  sampled  by  the  DPO  and  averaged  by  the  minicom- 
puter.  The  time-averaged  incident  waveform,  shown  in 
Figure  3  3a,  is  the  approximation  to  the  Gaussian  waveform 

described  by  (4-5).   The  signal  mean,  t    ,  is  two  nano- 
J  to  max' 

seconds,  and  the  standard  deviation,  a,  is  0.75  nanoseconds. 
To  measure  the  augmented  waveform,  the  target  was  placed 
directly  between  the  transmitting  and  receiving  antennas; 
a  distance  of  150  cm  from  the  transmitting  antenna,  and  a 
distance  of  12  0  cm  in  front  of  the  TEM  horn  antenna.   5  0 
augmented  waveforms  were  measured  and  averaged  exactly  as 
the  incident  waveforms  were.   The  time-averaged,  augmented 
waveform  is  shown  in  Figure  3  3b. 
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When  the  incident  waveform  was  subtracted  from  the 
augmented  waveform,  the  smoothed  impulse  response  for  the 
wire  scatterer  was  obtained.   This  measured  voltage  is 
presented  in  Figure  34,  and  is  directly  proportional  to  the 
scattered  electric  field  from  the  target.   Because  of  this 
relationship,  a  direct  comparison  between  the  measured 
voltage  of  the  target  response  and  the  theoretical  calcu- 
lations of  the  scattered  electric  field  can  be  made. 

The  method  of  moments  described  by  Harrington  in 

reference  16  is  the  basis  for  the  numerical  computation  of 

the  electric  field  scattered  by  the  wire.   In  the  computer 

program  used,  a  3.0  volt/meter  electric  field  with  Gaussian 

parameters  of  t    =2.0  nanosec  and  a=0.75  nanosec  was 

max 

broadside  incident  to  the  wire.   The  wire  was  divided  into 
22  segments  for  computational  purposes.   Scaling  the 
theoretical  values  obtained  by  a  factor  of  0.002  307  allowed 
the  theoretical  results  to  be  plotted  on  the  same  coordinate 
axes  as  the  experimentally  determined  values.   Agreement 
between  the  two  sets  of  data  is  seen  to  be  excellent  in 
Figure  34.   Of  the  13  theoretical  values  plotted,  only  one 
point  deviated  significantly  from  the  measured  data.   This 
occurred  at  the  first  negative  peak  of  the  response.   The 
deviation  arises  from  the  difference  between  the  actual  shape 
and  the  assumed  Gaussian  shape  of  the  incident  field  on  the 
wire. 
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The  close  correspondence  between  the  theoretical  and 
experimental  data  validates  the  measurements  made  in  the 
scattering  laboratory.   However,  further  calibration  of  the 
scattering  range  is  required  to  determine  the  amplitudes 
of  the  electromagnetic  fields  involved  in  scattering 
interactions.   This  will  be  necessary  to  develop  target 
classification  software  routines  in  follow-on  research. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  in  developing  the  time  domain  measurements 
scattering  laboratory  is  to  create  a  facility  which  will 
support  research  in  the  area  of  transient  electromagnetics. 
Ultimately,  this  will  assist  in  determining  the  feasibility 
of  the  design  of  advanced  radar  systems  based  on  principles 
of  time  domain  signal  analysis.   The  work  thus  far  com- 
pleted has  established  a  foundation  in  this  area  of  research 
and  has  provided  many  opportunities  for  continued  research 
in  greater  depth  and  on  a  broader  basis.   At  the  same  time 
new  problems  have  been  created  which  must  be  addressed. 

To  summarize  the  accomplishments  achieved  since  the 
project  conception;  there  is  now  available  a  complete 
physical  laboratory  which  contains  the  essential  electronic 
hardware,  antennas,  and  associated  equipment,  all  properly 
interfaced,  to  make  direct  time  domain  measurements  of  the 
scattered  electromagnetic  fields  from  mirror  imaged,  sym- 
metric objects.   There  exists  software  which  is  capable  of 
managing  the  signal  acquisition  process  of  the  electronic 
hardware,  reducing  the  noise  in  the  signal,  extracting 
the  smoothed  impulse  or  gaussian  pulse  response  of  a 
target  from  the  acquired  signal;  and  finally,  the  software 


91 


is  capable  of  storing  and  retrieving  any  necessary  data  or 
additional  program  instructions  for  further  signal  pro- 
cessing in  600  kilobytes  of  supplementary  tape  memory. 

A.   ANTENNAS 

A  problem  that  has  arisen  concerns  the  transmitting 
antenna.   Above  the  image  plane,  a  twenty-one  foot  wire, 
monopole  antenna  has  been  constructed  to  transmit  the  inci- 
dent electromagnetic  field.   Its  increased  length  over  a 
previous  antenna  tested  offers  a  longer  time  window  on  the 
oscilloscope  to  view  target  responses.   The  impedance  mis- 
match at  the  feedpoint  of  this  antenna  is  a  problem  worth 
further  investigation.   Although  the  reflections  at  this 
mismatch  are  now  within  the  system  tolerances,  this  was  not 
the  case  with  the  previous  antenna.   Reflections  from  its 
feedpoint  caused  occasional  pulse  generator  malfunctions. 
A  possible  recommended  approach  to  further  improve  the 
present  impedance  mismatch  of  the  monopole  is  to  build  a 
conical  outer  conductor  from  the  end  of  the  coaxial  cable 
to  the  image  plane  that  will  provide  a  gradual  transition 
in  the  impedance.   This  is  discussed  in  reference  12. 

In  addition  to  the  transmitting  antenna,  a  satisfactory 
receiving  TEM  horn  antenna  was  constructed.   Though  scatter- 
ing measurements  obtained  while  using  this  horn  antenna 
are  in  agreement  with  theoretical  results,  it  was  felt  in 
later  stages  of  the  experimental  work  that  the  TEM  horn 
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antenna  provided  more  signal  distortion  than  originally 
anticipated.   This  is  partially  confirmed  by  the  results 
of  the  impulse  response  of  the  coupled  antennas.   A 
thorough  experimental  analysis  is  recommended  along  with  a 
determination  of  how  to  optimize  signal  reception  on  the 
image  plane. 

B.   SIGNAL  PROCESSING 

Enhancing  the  signal  received  may  be  possible  with 
further  refinements  in  the  software  routine.   The  program 
in  Appendix  A  is  excellent  for  initial  signal  processing. 
This  includes  sampling  waveform  data,  incident  and  augmented 
signal  averaging,  and  determining  target  responses.   The 
graphical  display  presents  useful  waveform  comparisons 
along  with  correct  statistical  data.   Although  the  re- 
mainder of  the  program  has  been  proven  to  be  mathematically 
correct,  the  final  presentation  of  the  impulse  response 
of  a  target  is  polluted  with  so  much  noise  as  to  be  unusable 
in  most  cases.   The  causes  for  this  should  be  investigated. 
One  possible  cause  may  be  the  lack  of  use  of  functions  such 
as  the  Hamming  window  to  reduce  aliasing  in  the  FFT's.   The 
signal  processing  ROM  #2  user  manual  discusses  the  use  of 
the  "TAPER"  command  in  time  windowing  functions. 
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C.   TARGETS 

In  addition  to  the  targets  already  constructed,  several 
other  configurations  are  needed.   One  example  is  a  large 
hemisphere  made  from  aluminum.   This  will  serve  as  a 
calibration  tool  in  scattering  measurements.   Another 
example  is  a  set  of  bisected  cones  to  match  the  dimensions 
of  the  cones  already  fabricated.   These  bisected  cones  will 
provide  another  aspect  angle  and  incident  polarization  from 
which  to  view  scattering  from  this  particular  geometric 
configuration . 

In  general,  the  targets  should  be  larger  in  dimension. 
For  some  of  the  original  targets,  their  responses  are  less 
than  three  nanoseconds  duration  and  hence  are  indistin- 
guishable from  noise  in  the  system.   Sampling  by  the  DPO 
should  not  be  a  problem  related  to  target  resolution  because 
of  the  high  sampling  rates  used.   (25-30  picoseconds/sample) 
The  problem  is  dependent  more  on  the  amplitude  of  the  target 
response  than  its  duration.   Use  of  the  succeeding  pulse 
generator  may  serve  to  correct  this  problem  as  discussed 
in  Chapter  III.   Improvements  in  antenna  reception  and 
system  noise  reduction  should  also  contribute  to  improved 
target  resolution.   This  should  be  an  influencing  factor  in 
the  optimization  of  the  TEM  horn  antenna. 
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D.   CONCLUDING  REMARKS 

Future  plans  for  the  development  of  the  scattering 
laboratory  include  conversion  of  the  antennas  from  a  bistatic 
system  to  a  monostatic  system.   The  basic  monostatic  antenna 
design  is  patterned  from  information  found  in  reference  17. 
It  is  recommended  that  a  thorough  experimental  analysis  of 
the  transmitting  and  receiving  characteristics  of  this  conical 
antenna  be  made  prior  to  any  scattering  measurements  to 
insure  a  thorough  understanding  of  the  system.   It  is  also 
recommended,  in  general,  that  a  complete  knowledge  of  all 
aspects  of  the  scattering  range  be  understood  to  insure  that 
a  proper  interpretation  of  experimental  results  is  made  at 
all  times . 
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APPENDIX   A 
Software    Program 
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APPENDIX  B 
Target  Dielectric  Material  Specifications 

Properties  and  Applications  Data 

Dielectric  Constant,  60  Hz  to  12  GHz  2.54 

Loss  Tangent,  60  Hz  to  10  GHz  0.0005 

Homogeneity  of  Dielectric  Constant  +0.01 

Isotropy  of  Dielectric  Constant  +_0.01 

Volume  Resistivity,  ohm-om  >1016 

Surface  Resistivity,  ohms/square  1014 

Dielectric  Strength,  volts/mil(Kv/mm)  500(20) 

Specific  Gravity  ~  1.06 

Tensile  Strength,  psi(Kg/cm  1  11,000(770) 

Flexural  Strength,  psi(Kg/cm  )   2  17,000(1,190) 

Modulus  of  Elasticity,  psi(Kg/cm  )  300,000(21,000) 
Izod  impact,  ft.  lbs. /in.  of  notch 

(Kg-cm/cm)  0.4(2.2) 
Coefficient  of  Linear  Expansion/°C  fi       R 

(/°F)  63x10   (35x10   ) 

Thermal  Conductivity,  ( BTU) ( in) / (hr )  0.87 

(ft2)(°F)(cal)(cm)/(sec)(cm2)(°C)  0.0003 

Rockwell  Hardness,  M  Scale  105 
Water  Absorption,  %  gain  in  24  hrs . 

at  25°C  0.003 
Recommended  maximum  operating 

temperature : 

Continuous  125°C(257°F) 

Short  Time  200°C(392°F) 
Recommended  minimum  operating 

temperature  -70°C(-94°F) 
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